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Abstract 

We present a critical analysis of Pade-based methods for the unitarization of low energy 
amplitudes. We show that the use of certain Pade Approximants to describe the resonance 
region may lead to inaccurate determinations. In particular, we find that in the Linear Sigma 
Model the unitarization of the low energy amplitude through the inverse amplitude method 
produces essentially incorrect results for the mass and width of the sigma. Alternative sequences 
of Pades are studied and we find that the diagonal sequences (i.e., [N/N]) have much better 
convergence properties. 



1 Introduction 

Effective Field Theories (EFT) have become a very useful tool for the description of low-energy 
physics PP. Based on symmetry and dimensional analysis arguments, they allow to organize a 
perturbative expansion of the amplitudes in powers of soft momenta over some characteristic scale, 
p/A. For instance, below the lightest resonance multiplet, the interactions between the pseudo- 
Nambu-Goldstone bosons from the spontaneous chiral symmetry breaking -the pions in the SU(2) 
case- is provided by Chiral Perturbation Theory (xPT) [21 El S] - Its range of validity is a few 
hundred MeV, breaking down as one approaches new states not included in the EFT. Thus, xPT 
is unable to describe the resonance region and one needs to incorporate additional ingredients. 

Several works have tried to extend the range of validity of the EFT by means of the unitarization 
of the low-energy amplitude. Unitarization methods have been applied extensively to Quantum 
Chromodynamics and xPT [H]> where the issue of the scalar resonances is particularly interesting. 
However, their range of applicability is wider; for example they have also been applied to WW- 
scattering [6j. Pade Approximants (PAs) [8] and the Inverse Amplitude Method (IAM) [9l 110] 
are among the most usual ones, although some remarks and criticisms on the reliability of these and 
other unitarization methods have been raised [U \TT\ [T2] . In the present letter we focus our attention 



on the PAs. We study up to what point one can rely on them to describe the resonant properties of 
the theory. By means of a couple of models (the Linear Sigma Model and a vector resonance model) , 
we show that these militarization procedures may lead to improper determinations of the resonance 
pole position (masses and widths) [13] . Furthermore, one may not recover the right values for the 
Low-Energy Couplings (LECs) of the EFT if the PAs are applied to describe the resonant region. 

The starting point of our analysis is, therefore, a model where the properties of the resonances 
are known. Then we derive and unitarize the corresponding low-energy amplitude. The predictions 
for the masses and widths obtained from the PA sequence [1/iV] (referred to as IAM in some 
works [8| \9\ \10\ I12j) are compared and found to be quite different from those of the original model. 
As an alternative to the badly behaved sequence [1/iV], we propose the use of the PA sequence 
[N/N], which quickly converges as iV is increased. For simplicity the chiral limit will be assumed 
all along the paper, but this will not alter the main conclusion. 

The structure of this letter is the following. We begin, in Section [21 with the analysis in the LSM: 
in Section [2. II we compute the exact position of the pole of the sigma correlator at the one-loop level. 
Then in Section [2.21 we consider the Linear Sigma Model (LSM) at low energies and present the tttt 
scattering amplitudes. In Section [2731 we apply the IAM procedure to the partial wave amplitudes 
to recover the mass and width of the sigma. In Section [3] we analyze PA sequences at higher orders. 
In Section 2] the corresponding results are discussed in the context of a vector resonance model. We 
then conclude, in Section \5\ with a short discussion. 



2 One-loop Linear Sigma Model 



2.1 Sigma pole position up to O(g) 

We begin presenting the results for the one-loop corrections to the mass and the width of the 
sigma in the LSM. At tree-level, the sigma mass is found to be M 2 = 2// 2 , and the width is zero. At 
next-to- leading order, the sigma pole gets shifted due to the quartic potential, i.e., M 2 = 2/u 2 + 0(g), 
and the width becomes different from zero. 

In order to determine the scalar meson mass and width up to O(g), we compute the one- loop 
sigma correlator [3], 
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where p(s) = \J\ — AM 2 /s and the term —13/3 is determined by the renormalization scheme 
chosen by Ref. [3] , which sets the relation 2gF 2 = M 2 at the one-loop order, with F the pion decay 
constant. Now it is possible to extract the pole s p of the propagator up to the considered order in 
perturbation theory. If one approaches the branch cut from the upper part of the complex s-plane, 
the pole in the second Riemann sheet is located at 
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where we have used s p = M 2 + 0{g). The pole mass and width, defined from s p = (M p 
are then given by 
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As expected, at lowest order the pole width agrees with that derived from the decay amplitude 

2.2 Low-energy expansion 

We now consider the LSM at low energies. The contribution from the sigma exchanges to the 
renormalized 0(p 4 ) couplings gives [3] 
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The 7T7r-scattering is determined by the tt + tt — > irir amplitude, which is given up to 0(p ) in 
the chiral expansion by 
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where /i refers here to the arbitrary renormalization scale, and the chiral limit has been considered. 
With this one constructs the definite isospin amplitudes 
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The partial wave projection is then provided by 

tj{s) = -^J dco&O Pj(cos9) T{s,t,u)\ 
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where 9 is the scattering angle in the irir center-of-mass rest frame. 

Hence, for the first partial waves with IJ = 00,11,20, one finds the following 0(p 2 ) 

amplitudes, 
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and at 0(p 4 
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For the next section, it will be suitable to rewrite the 0(p 4 ) amplitudes in terms of the LSM 
parameters and the 0(p 2 ) scattering: 
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where we have used Eq. and the relation 2gF 2 = M 2 [3] . 

2.3 Unitarization of the %PT amplitude 

The Inverse Amplitude Method (IAM) provides an amplitude that is unitary not only at the 
perturbative level but exactly. This means that in the elastic limit one has for s > the partial- 
wave relation 

lmt(s) = \t(s)\ 2 , (11) 

where the indices I J are assumed (t = tj). This relation can be reexpressed as a relation for the 
inverse amplitude: 

Imt(s)" 1 = -1. (12) 

Thus, the imaginary part of i(s) -1 becomes completely determined and one only needs to spec- 
ify the real part Ret" 1 . The IAM relies then on a low-energy matching to xPT (with t~£. 
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1 — i(4)A(2) + ••• ) m order to fix the unknown part of the amplitude. Thus, at 0(p 4 ), one has 
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the unitarized amplitude, 
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This expression is sometimes also known as a P 1 Pade Approximant of the partial-wave amplitude. 
The IAM has been also extended up to 0(p 6 ) by means of what is sometimes named as a P 2 l 
approximant [10l 02] : 
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However, we want to remark that t(s) IAM is not a PA in the variable s: It is not a rational 
approximant since it also contains the logarithms from the pion loops. Thus, strictly speaking no 
theoretical argument ensures the recovery of the physical amplitude. Only in the tree-level limit 
t(s) IAM becomes a PA. In any case, we will see that both the whole and the tree-level IAM amplitudes 
are unable to reproduce the original partial waves in the resonance region. 

Given the 0{p 2 ) and 0(p A ) %PT amplitudes from Eqs. (i8l)- (fT0l) . it is then possible to extract 
the poles of the corresponding t(s) IAM for the LSM, satisfying 1 = t(s)u\/t(s)r 2 ) at s = s p : 
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These are the poles that appear in the unphysical Riemann sheet as one approaches from upper half 
of the first Riemann sheet. There is also a conjugate pole at s* if one approaches the real s-axis 
from below. 

The first thing to be noticed is that poles appear in the I, J =11 and 20 channels even for small 
values of g, contrary to what one expects in the LSM, where no meson with these quantum numbers 
exists. Furthermore, these "states" are not resonances, as they are located on the left-hand side of 
the complex s-plane, out of the physical Riemann sheet, and carrying a negative squared mass. 

As for the IJ = 00 channel, one finds a resonance with pole mass and width, 
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The IAM predictions for M 2 and M p T p result, respectively, 40% and 80% smaller than the original 
ones in the LSM, computed in Eq. ([3]). This points out the low reliability of this particular method 
in order to recover the hadronic properties of the theory from its effective low-energy description. 

The IAM poles remain badly located even in the weakly interacting limit, so this failure cannot 
be attributed to non-perturbative effects. In the limit when g — > and M a is kept fixed one finds 
that the poles predicted in all the different channels fall down to the real s-axis. We are left with 
just tree-level amplitudes and the expressions become greatly simplified. Due to the smoothness of 
this limit, it will be assumed in the next analysis of higher order Pade Approximants [M/N] and 
in the study of the vector model in Section HI 
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3 Higher order Pade Approximants for tree-level amplitudes 



In this section we consider higher order Pade Approximants to the partial wave amplitudes, with 
the hope that this will provide some insight on the nature of the unitarization process discussed 
above. We will see that the PA sequence associated with the IAM does not converge properly, and 
that diagonal sequences are much more suitable for this purpose. We begin with a brief overview 
of the theory of PAs. 



3.1 Generalities on Pade Theory 

Let /(s) be a function with a Taylor expansion around s = 0. A Pade Approximant (PA) to f(s), 
denoted as Pff (s), is defined as the ratio Ppf (s) = Rm(s)/Qn(s) of two polynomials of orders M 
and N, respectively, and such that f(s) — Pjf(s) = 0{s M+N+1 ) when s — > 0. 

The convergence properties of the PAs to a given function are more complex than those of 
the Taylor expansion. However, they converge for a broader set of functions, even in the case of 
slowly convergent or asymptotic power series, and they usually carry smaller errors than the Taylor 
expansions (when these are applicable). Furthermore, in many cases, the Pades have been found to 
provide a fairly good approximation even beyond their expected range of applicability. 

Pommerenke's theorem states that the sequence of diagonal PAs, i.e. [N/N], to a meromorphic 
function is convergent everywhere in any compact set of the complex plane except, at most, for a zero 
capacity set PS E51 E3 ESI E31 EBJ. This obviously includes the poles of /(a), where the original 
function is ill-defined. In addition, the PA may produce a series of poles absent in f(s). Thus, for a 
given compact region K, in the complex plane, Pommerenke's theorem of convergence requires that, 
either these undesired poles move away from the region fC as the order of the PA increases, or they 
pair up with a close- by zero becoming what is usually called a defect [1 16] . Although the PA breaks 
down in the very neighborhood of these extraneous poles, away from them the approximation is 
safe. Likewise, a PA can approximate a multivaluated function (for example, a function with a 
logarithmic branch cut). The sequence of PAs will cluster its poles along the cut [16j . as we will see 
in the following section. 



3.2 Tree-level PAs in the LSM 



In order to be able to handle the amplitude at higher orders, we will consider the tttt scattering 
at tree-level. This is equivalent to working in the limit g < 1 and keeping just the first non-trivial 
contribution in the g expansion. Thus, the 7T7r-scattering is determined in the LSM by the function 
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In addition to meromorphic, there are known theorems of convergence with PAs for Stieltjes functions, continued 
functions, Gauss hypergeometric functions, Bessel functions, some kind of divergent series, sets of complex points, 
etc. This has been applied in the past to various kinds of scattering processes |16| . 
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The In [l + s/M 2 ] logarithms come from the partial- wave projection of the tree- level exchanges of 
resonances in the crossed channel. They have absolutely nothing to do with the logarithms of the 
xPT amplitudes in Eq. (jlOp . which come from the tttt loops. 
At low energies the amplitude becomes 
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so the partial waves are given by, 
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The comparison between the low-energy expansions and the whole result provides a first insight of 
the piece of information that is lost in the unitarization procedure. At high energies, the partial 
waves contain poles on the right-hand side of the s-plane, related to s-channel resonance exchanges, 
and a left-hand cut, related to the crossed-channel resonance exchanges. At low energies, both kinds 
of exchanges contribute equally to the low energy couplings, so the crossed resonance exchanges shift 
the IAM poles from their physical value. Although t and u channels are not so relevant in the region 
close to the resonance pole, at low energies they are as important as the s-channel. 
The simplest Pade, P±, gives the prediction 
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which agrees with the one-loop calculation from Eqs. (|14j> — (|16 p if one remains at leading order in g. 



3.3 Higher order PAs in the LSM 

The convergence of a sequence [M/N] of Pade Approximants to a function implies that the PA 
tends to mimic its analytical structure as M,N — > oo. The common procedure for the construction 
of a sequence of PAs is to increase M, N following a given pattern, e.g. Af = N — > oo. In some cases, 
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Figure 1: Left: Position of the nearest pole to M% for the first PAs of the form [1/N] with ./V odd (for even 
N all the poles are complex). Right: Poles of the Pg X in the complex plane. 



this allows the application of known mathematical theorems that ensure convergence |15[ [T6| [20] . 
This has been used for the study of certain Green functions |17[ [TBI H9j - However, little is known 
about the sequence [M/N] when M is kept fixed and N — > oo (for instance Pf , P\ ...). 

We will employ and compare the [1/N] and [N/N] sequences for the study of the tttt partial 
wave scattering amplitudes. We will also comment on PAs of the [N + K/N], e.g [N — 2/N]. In 
the next lines we will focus our attention on the 13 = 00 partial wave, but analogous results are 
found for the other channels. Former works pointed out that the PAs and other unitarizations fail 
to incorporate the crossed channel resonance exchanges [141 124] . Nonetheless, we will see that as 
N grows, the poles of the sequence [N/N] actually tend to mimic not only the s-channel poles but 
also the left-hand cut contribution from diagrams with resonances in the t and u channels. 

The sequence [K/N], with K fixed, is studied in the present section in the framework of the 
Linear Sigma Model. It provides an example of the behavior of these kind of sequences. K = 1 is 
chosen because of the similarity of this sequence and the IAM [9j [T0l [T2l [131 124] . Our results are 
summarized in Fig. [TJ No convergence is found with this sequence. In the case of N odd, Fig. [TJa. 
shows that the Pjy pole closest to M 2 does not approach this value even for very large N, always 
remaining a 30% below. The analytical structure of the original amplitude (s-channel sigma pole 
plus left-hand cut) is never recovered since the [1 /N] PAs always set the poles in the circular pattern 
shown in Fig. [TJb. This suggests that the use of further [1/N] approximants to extend the IAM is 
not the optimal way to proceed, even if we had an accurate knowledge of the low-energy expansion 
up to very high orders. 

Alternatively, the use of sequences such as [N + K/N] (e.g. [N-2/N], [N — 1/N], [N/N], [N + 
1/N]...) seems to be a better strategy. In the following we analyze the sequence [N/N], as it 
ensures the appropriate behavior at high energies, \t(s)\ < 1. Nevertheless, similar results have 
been generally found for the [N + K/N] PAs with if The Pjy pole closest to M% is shown in 
Fig. [2ja. One finds a quick convergence of the sequence: P* reproduces the sigma pole a 40% off but 
P 2 2 disagrees by less than 1%, Pf by less than 0.1%, etc. Notice that already P| provides a much 
better description than Pq\, although one includes far more low-energy information in the latter. 
All this points out the sizable discrepancy of the first element of the sequence (P/) with respect to 
the original amplitude. It also indicates that the [1/N] PAs do not produces a serious improvement. 
On the contrary, the [N/N] sequence provide a far more efficient strategy with a quick convergence. 

Likewise, Fig. [TJb shows how the [1/N] PAs are unable to recover the analytical structure of the 
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original amplitude, whereas the [N/N] sequence, besides providing the isolated pole of the sigma, 
tends to reproduce the left-hand cut as N increases. The poles of P|o are plotted in Fig. [2 Although 
a PA is a rational function without cuts, these are mimicked by placing poles where the cuts should 
lie. The Pfo nas one isolated pole near Ml (with an accuracy of 10 30 ) and nineteen poles over 
the real axis at s p < —Ml, i.e. on the left-hand cut of the original function. As N is increased, the 
number of poles lying on the branch cut increases too. 

A remarkable feature found for the first P$ approximants (P*, Pf , P3) is that they obey exact 
unitarity, as it happened with the IAM sequence [1/iV]. 

We would like to remark that, although we lack of a general theorem that ensures the convergence 
of the [N + K/N] sequences, we have found that they reproduce the original partial waves for 
arbitrary K. Moreover, after performing modifications on the structure of the amplitudes of Eq. (|19[) 
we still found convergence the different [N + K/N] PAs. In several situations the PAs set all poles 
over the left-hand cut position and one isolated pole that approached Ml when N — > 00. In the 
worst cases, in addition to this we found extraneous poles that either moved away as A^ increased 
or they tend to be canceled by nearby zeros at A^ —* 00. 

As an amusement, we have also probed the PA sequence [N/1] which has the same number of 
inputs as the [1/N] has for a given N. In this new case we have found convergence in both LSM and 
the resonance model presented in the following section but slower than the [N/N]. For instance, 
the prediction for the Ml for the first P^ are 4%$ = 0.55, 1.47, 0.73, 1.27, 0.81... A criticism that can 
be done to this sequence is its lack of unitarity, in contrast to the other studied sequences. 



4 Vector Resonance Model 



In order to broaden our analysis, we consider now a model with just vector mesons |14| . It could 
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be derived either from the gauged chiral model [2T] for the couplings 3g p F 2 = M 2 , or from resonance 
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The study of the I J =11 partial wave leads to the same conclusions found for the LSM. It is 
at first sight remarkable that, on the contrary to the previous case, one recovers s p 



M 2 from the 
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first-order approximant P/. However, the sequence [1/N] already worsens at N = 2, where the two 
complex-conjugate poles are located at s p = (0.71 ± 0.96i) Mp 2 on the physical Riemann sheet. On 
the other hand, P$ exactly recovers s p = Mp 2 for any odd N. For N even, the prediction from P| 
is a 30% off, but one has again a quick convergence to s p = Mp 2 as N increases: Pf disagrees by 
less than 0.1%, Pg disagrees by less than 10~ 6 , etc. 

Furthermore, the [1/N] and [N/N] sequences produce, respectively, the same structure of poles 
found for the LSM. This is, [1/N] generates the circular structure of poles of Fig. [TJb. and the 
sequence [N/N] places one pole at s p ~ Mp and the remaining ones reproducing the left-hand cut 
in analogy to Fig. [2jb. 

5 Discussion 

In this letter, we have addressed the reliability of the unitarization of low-energy amplitudes 
through Pade Approximants. The one-loop analysis of the LSM has led to IAM predictions of the 
sigma mass and width respectively a factor 2 and 5 smaller than the original ones in the model. 

Some disagreement has been found in phenomenological determinations of the /o(600) pole 
width |llj . pointing out that the IAM is unable to recover at the same time the value for the mass 
and the width predicted from Roy equations |23j. This discrepancy becomes even more obvious in 
the chiral limit, with differences much larger than the expected O(m^) corrections. 

It has been also shown that the P* mass prediction from the one-loop analysis agrees at leading 
order with that coming from the tree-level amplitude. This argument allowed the construction 
of higher order PAs based on the tree-level amplitude of the LSM and a vector model. It was 
found that the sequence [1/N] is unable to reproduce the original partial wave scattering amplitude 
whereas the PAs of the form [N/N] display a quickly convergent behavior. Moreover, though we 
lack of convergence theorem, it was found that the [N + K/N] PAs were able to reproduce the 
partial waves for arbitrary K. 

For all this, we suggest the use of the [N/N] sequence rather than [1/iV]. Unfortunately, the 
study on broad resonances requires to go beyond the tree-level approximation. Thus, in the real 
world, the I J = 00 channel needs the inclusion of loops, precluding by now the extension to PAs 
beyond 0(p 4 ) in the chiral expansion, i.e. P*. There is also a clear limitation on our experimental 
knowledge of the low-energy couplings, which barely goes beyond 0(p 4 ). However, our proposal 
should be still suitable for the analysis of theories with narrow resonances and a relatively good 
knowledge of the experimental low-energy amplitudes. 
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